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ABSTRACT 

The properties of galaxies in the local universe have been shown to depend upon their 
environment. Future large scale photometric surveys such as DBS and Euclid will be vital to 
gain insight into the evolution of galaxy properties and the role of environment. Large sam¬ 
ples come at the cost of redshift precision and this affects the measurement of environment. 
We study this by measuring environments using SDSS spectroscopic and photometric red- 
shifts and also simulated photometric redshifts with a range of uncertainties. We consider the 

nearest neighbour and fixed aperture methods and evaluate the impact of the aperture pa¬ 
rameters and the redshift uncertainty. We find that photometric environments have a smaller 
dynamic range than spectroscopic measurements because uncertain redshifts scatter galaxies 
from dense environments into less dense environments. At the expected redshift uncertainty 
of DBS, 0.1, there is Spearman rank correlation coefficient of 0.4 between the measurements 
using the optimal parameters. We examine the galaxy red fraction as a function of mass and 
environment using photometric redshifts and find that the bivariate dependence is still present 
in the SDSS photometric measurements. We show that photometric samples with a redshift 
uncertainty of 0.1 must be approximately 6 — 16 times larger than spectroscopic samples to 
detect environment correlations with equivalent fractional errors. 

Key words: galaxies: evolution - galaxies: formation - galaxies: clusters: general - galaxies: 
photometry - galaxies: statistics 


1 INTRODUCTION 


Fluctuations in the early universe were the seeds of the structur e 
that is observed today. After a period of rapid inflation (lGuthll98ll) . 
in some regions, under the influence of gravity, d ark matter ac¬ 
cumu lated and subsequently collapsed into haloes dWhite & ReesI 
Il978h . This in turn, deepened the gravitation potential in these re¬ 
gions. In the hierarchical model of strueture growth; haloes grew 
by accreting dark matter and also through merging with satellite 
haloes. Baryons fell into the potential wells of the haloes and ev en- 
tually galaxies formed and evolved dPress & Schechtejll974h . It 
seems plausible that there should be a measurable connection be¬ 
tween the properties of galaxies and their environments. 

It is for this reason that galaxy environment has been an ac¬ 
tive research fleld for several decades. The pioneering work in 
this fleld was the discovery of the galaxy morphology-density re¬ 
lationship (IOemlej[l974 IPresslej 198Clf) . This work showed that 
elliptical galaxies are preferentially found in dense environments, 
whereas spiral galaxies are preferentially found in less dense envi¬ 
ronments. Large vol ume galaxy surv eys such as the Sloan Digital 
Sky Survey (SDSS) (lYork et al .120001) . utilizing modern CCDs have 
revolutionised the study of galaxy environment. Nearly complete 
galaxy catalogues, including the angular and redshift information 
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for large fractions of the sky, for the local universe, have made it 


possible to study the environmental dependence of simple galaxy 

properties, such as colour jHogg et all 12004: 

Blanton et al.ll2005l: 

Chuter et al. 1201 lb, morphology. JPark et al. 

I 2 OO 7 I; iKovac et alJ 

2010l:lBamford et alJl2009l: Skibba et alJl2009b 

1 star formation rates 

1 stellar populations 
Ltistieally signiflcant 

(iKauffmann et al.l 2004: Peng et al. 

2010b anc 

dThomas et al. 2010l; iRogers et al.l 

201(]|). Sta 


samples can be obtained from today‘s datasets, even after the 
data has been binned on several axes. A ga l axy bimodality has 
been firmly established dStrateva et al] l200ll : iBaldrv et^ |2004 
iMignoli ^ alJ2009l:IPozzetti et alllOlOh . There is a “red sequence” 
population of galaxies which tend to be redder, more massive, ellip¬ 
tical and live in more dense environments. Th e stellar populations 
in these galaxies formed on shorter timescales (iThomas et al .l2005h 
and the stars evolved rapidly until the gas and dust was depleted. 
After this their evolution was passive. Conversely the “blue cloud” 
population of galaxies tend to be bluer in colour, less massive, have 
spiral morphologies and live in less dense environments. They have 
evolved more slowly and still have active star formation. 

Both mass and galaxy environment are considered as drivers 
of galaxy evolution. However their precise roles and relative im¬ 
portance have been strongly debated. The deflnition of the mass 
and environment is vital. Galaxy environment is considered to be a 
sec ond order driver but i s potentially important for low mass galax - 
ies (iThomas et al.ll201(il : [Rogers et al.ll201Ql : IVulcani et afll^ll) . 
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The fraction of red galaxies, in samples from the SDSS, has been 
shown t o be a bivariate function of stellar mass and projected local 
density (iBaldrv & Baloghl[2006l : IPeng et al.ll201C)h . Hydrodynam- 
ical and semi analytical simulations, however, suggest that halo 
mass is the fundamental parameter. Observationally it is difficult 
to disent angle the effects of halo mass from galaxy environment, 
although baas et all (12012h have supplied a possible approach. Ev¬ 
idence that galaxy environment plays a signi ficant role in driving 
the evolution of satellite galaxies is mounting (van den Bosch et all 


I 2 OO 8 I: lFaltenbachei]l201cl:IPeng et all2012l: ICarollo et alJl2013h~ 

Finding correlations between the properties of galaxies and 
their environments does not, however, imply causation. It is dif¬ 
ficult to isolate and quantify, in a statistically robust manner, the 
relati ve importance of physical me chanisms such as tidal s trip- 
ping dRead et al.ll2006h. harassmen t dFarouki & ShaDirolll98lh and 
strangulation dLarson et al.lll98C)h that transform galaxy proper¬ 
ties, driven by a galaxy‘s environment. Ho wever there is substan¬ 
tial evidence for ram-pressure stripping d Abramson et al.l I 2 OIII : 
iBoselli & Gavazzil I20l4 iFumagalli et al.l 1^1 4) in local clusters 
obtained from multi wavelength observations. Dwarf galaxies in 
falling i nto clusters can be completely stripped of their gas in 
~1 Gyr dKennev et 

To study the redshift evolution of galaxy properties as a func¬ 
tion of environment statistically significant samples of galaxies are 
required for a large redshift range. The current state of the art 
surve ys are pencil beam surveys such as pEEP2 (Newman etalJ 
l2013h . UKIDSS Ultra deep survey tUPy dChuter et alJl201ll) and 
the COSMOS survey llKoekemoer et alJl2007h . These surveys are 
complete to relatively high reds hifts so facilitate stu dies of the evo¬ 
lution of galaxy properties te.g. lCooDer et al.|[2007h . However they 
target galaxies within relatively narrow areas (a few squared de¬ 
grees). Larger areas (thousands of square degrees) are required to 
establish galaxy samples that can be investigated further (e.g. by 
binning in mass, galaxy type, environment etc.). This is feasible 
only with photometric surveys. Currently, spectroscopic surveys 
with these requirements would be too costly and slow. 

The next generation of large scale photometric surveys in¬ 
cludes the Dark Energy Survey (DES) and ESA‘s space mission 
Euclid. The DES aims to survey 5000 square degrees of the sky 
by imaging 300 million galaxies up to a redshift of ~1.5 whereas 
Euclid will be a full sky survey to a redshift of ~2.0. These sur¬ 
veys will provide the opportunity to study the evolution of galaxy 
properties as a function of environment. Photometric surveys make 
it possible to survey an enormous number of galaxies as the av¬ 
erage measurement time per galaxy is considerably less than it is 
for spectroscopic measurements. This comes at the cost of redshift 
precision. In the DES the distribution of the difference between the 
photometric redshifts and the true redshifts is expected to approx- 
imate ly follow a Gaussia n distribution with a standard deviation 
~0.1 dBanerii et alJl2008h . 

Errors in the redshifts of the galaxies result in errors in the en¬ 
vironment measurements. A large error on a single redshift could 
result in a galaxy cluster member being measured in a void; a com¬ 
pletely different environment. The situation becomes worse when 
multiple galaxies in the same spatial vicinity are assigned redshift 
errors. The environment of a particular target can vary dramatically. 
Redshift information derived from photometric techniques will be 
available in the next generation of photometric galaxy surveys but 
spectroscopic measurements will be limited. The purpose of this 
paper is to investigate two commonly used methods to measure 
galaxy environment and explore the impact of using photometric 
redshifts. Previous work examining this theme with considerations 


of the DEEP2 survey includes ICooper et ^ (l2005h . Here we reex¬ 
amine the impact of photometric redshifts on galaxy environment 
measurements using SDSS data and consider the next generation 
of large scale photometric surveys. 

The paper consist of six sections. In the next section a brief 
review of the methods to measure galaxy environment is presented. 
Section [ 3 ] describes the data selection and galaxy properties used 
in this work. The results are described in Section|4|and further dis¬ 
cussion is presented in Section [5] We present the conclusions in 
Section [6| 

In this work we have assumed a cosmology with Vtrn = 0.3, 
= 0.7 and iTo = 70 kms“^Mpc“^. We adopt comoving coor¬ 
dinates to calculate the distances between galaxies as this system is 
suitable for this low redshift work {z < 0.1). 


2 REVIEW OF ENVIRONMENT METHODS 

In essence, a galaxy‘s environment is the density field in which a 
galaxy resides. We illustrate this in Figure [T] using SDSS data. The 
left hand plot shows the SDSS footprint and the right hand plot 
shows the average projected galaxy environment for the redshift 
range: 0.02 — 0.085. Several research fields, although often con¬ 
sidered as separate, probe this density field. For example, there is 
significant effo rt in the weak tensing co mmunity to produce den¬ 
sity maps (^e.g. ISzepietowski et alJl2014l) and the 2-point correla¬ 
tion function of the galaxy dist ribution is routinel y used in stud¬ 
ies of large-scale structure (^e.g. lToieiro et alJ[2014h . However, the 
environmental dependence of galaxy properties has been typically 
studied with two broad families of methods: the nearest neigh¬ 
bour methods and the fixed aperture methods. Within these two 
broad families of me thods there are a large number of variants (e.g. 
iMuldrew et al.ll201^ and this makes it challenging for researchers 
to compare results. The environment measurements computed via 
different methods should be correlated. However, they often probe 
different scales. For example, nearest neighbour methods in¬ 
trinsically adapt the scale probed for each target whereas the scale 
probed by aperture methods is fixed. 


2.1 Fixed aperture 


Perhaps the simplest method to measure galaxy environment is to 
count the number of galaxies (excluding the target) within a fixed 
volume centred on the target. A number density can be computed 
by dividing the number of galaxies found by the comoving volume. 
This type of approach is known as a fixed aperture me thod. A vari¬ 


ety o f apertures have been used inclu ding: spheres dCroton et al. 


2005 ). cylinders (iGallazzi et al.l |200^ and annuli (IWilman et al 


20101) . These fixed aperture methods are discrete measures of envi¬ 
ronment as the number count is either a positive integer or zero and 
the volume is fixed. The discrete nature of these environment mea¬ 
sures makes it difficult to distinguish between low density environ¬ 
ments. Field galaxies can exist in isolation with no other galaxies 
present within an aperture, centred on the galaxy. 

The size of the aperture used in a study, for example, the ra¬ 
dius of the sphere or the radius and height of the cylinder, must 
be chosen. Choosing a small aperture will result in a large num¬ 
ber of the target volumes being devoid of galaxies. However the 
measure would be sensitive to the most dense environments. Con¬ 
versely, a large aperture will reduce the number of target volumes 
that are devoid of galaxies, enabling measurements of a range of 
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Figure 1. Galaxy distribution for the SDSS main footprint (left) and the projected galaxy environment (right) for the redshift range: 0.02 < 2 ; < 0.085. Dense 
environments are shown in red and sparse environments are in blue. 


low density environments. Probing increasingly large scales, how¬ 
ever, tends to homogen ise the environinents. A large range of scales 
from <1 — 10 Mpc (iMuldrew et al.ll201^ have been studied in 
previous works probing individual haloes through to structures in 
the cosmic web. Studies show that environmental effects on galaxy 
properties (e.g. colour) are strongest on scales s maller than ~1 Mpc 
telanton & Berlindll2007l : Iwiman et al.ll2010l) . Ideally the choice 
of scale should match the scale of the physical processes that are 
believed to drive the evolution of the galaxy properties. However 
considerations of the sample size, signal-to-noise and the dynamic 
range of environments often leads to a compro mise, particularly 
when volume limited samples are constructed (iKauffmann et al.l 
l2004h . 


2.2 Nearest neighbour 

The other main family of methods is the nearest neighbour 
methods. For each target galaxy in the dataset the nearest 
neighbour g alaxy is identified. Spherical volumes can be used with 
simulations (iHaas et alJl2012l) as the distances between galaxies 
can be determined accurately. More often in observational studies 
cylindrical volumes are used to constrain the nearest neighbours, 
to mitigate the effects of redshift distortions or redshift errors. In 
either case the comoving distance to the nearest neighbour is 
computed and this is used as the radius to compute a comoving vol¬ 
ume centred on the target galaxy. Dividing N by this volume gives 
a number density. A dense environment is obtained when the A^* 
nearest neighbour is close to the target. 

A consequence of the A^* nearest neighbour method is that 
if there are A^ — 1 neighbours in a group, near the target, the A^* 
neighbour will be located in the next, closest group. Hence the vol¬ 
ume targeted will be based on the inter and not the intra group dis¬ 
tance. A careful choice of the value of N enables studies to focus 
on a particular range of scales. 


3 DATA SELECTION 

To examine the loca l unive rse we follo wed the app r oach adopted 
by iBaldrv & BaloghI (I 2 OO 6 I) and also bv IPeng et al.l (|201Q|) . Spec¬ 
troscopic and photo metric datasets were ex tracted from the seventh 
data release (DRV) (lAbazaiian et alJ[200^ of the SDSS using the 


CasJobs website. Specifically this data includes the SDSS spectro¬ 
scopic and photometric redshift measurements. 

First the photometric data was obtained from the “Galaxy 
View” in the SDSS database. Only objects with a Petrosian, r, mag¬ 
nitude in the range: 10.0 < r < 18.0 were selected and to ensure 
stars were exclude d only those objects where rpsF - r’Modei> 0.25 
were included (e.g. lScranton et alJ2002h . We refer to this sample as 
the photometric parent sample as all of the objects in the subsequent 
samples in this work can be found in this sample. Spectroscopic 
data was extracted for all of the objects in the photometric parent 
sample where measurements were available. The photometric par¬ 
ent sample contained 1,607, 820 objects of which 803, 939 objects 
had spectroscopic measurements. These selections were reduced to 
obtain two volume limited samples (see Section [TM . One sam¬ 
ple was based on the spectroscopic redshifts and one was based on 
the photometric redshifts. During this reduction the redshift range 
was restricted to: 0.02 < z < 0.085 and galaxies where the stel¬ 
lar mass and absolute magnitude were unavailable were excluded. 
This ensures that the photometric and spectroscopic samples were 
reasonably well matched. 

To investigate a range of redshift uncertainties, in addition to 
the spectroscopic and photometric data extracted from the SDSS, 
simulated photometric redshift catalogues were created (see Sec¬ 
tion [3]3]). Volume limited samples for the simulated redshift cata¬ 
logues were created using the same methodology (see Section [T4l ). 

3.1 Stellar masses, absolute magnitudes and k-corrections 

Stellar masses are adopted from (iMaraston et aLl201^ . These have 
been derived through a standard SED fitting procedure from 
the SDSS obser ved u,g,r,i and z band m agnitudes using the fitting 
code HYPERZ (iBolzonella et all[200^ based on the SDSS spec¬ 
troscopic redshifts. In this work we used th e output based on the 
the passive templates from iMarastonI (l2005l) consis ting of a metal- 
poor a nd a metal-rich old population as described in iMaraston et alJ 
(I 2 OO 9 I) . These are called “passive” because they exclude very recent 
and current star formation, but still allow for a wi de range of stel¬ 
lar pop ulation ages of 3 — 15 Gyr. It is shown in IMaraston et alJ 
(I 2 OI 3 I) that these are adequate to properly model the stellar pop¬ 
ulations also of (star forming) galaxies in the blue cloud as they 
minimise the effect of underestimating th e age (and hence th e 
galaxy mass) owing to recent star formation (IMaraston et alJ201^ . 
Hence these provide the most accurate mass estimates also for star 
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Figure 2. Sample selection in the Mr-redshift plane for the spectroscopic 
catalogue. The vertical magenta lines mark the redshift bounds of the selec¬ 
tion. The connecting magenta curve shows the absolute magnitude detec¬ 
tion limit. The horizontal magenta line binds the density dehning popula¬ 
tion. The number of galaxies in each 2-dimensional bin is represented with 
a colour as indicated in the colour bar. 

forming galaxies. We have verified that our masses are consistent 
with the publicly available MPA- jhlQ SDSS DR7 mass estimates 
(iKauffmann et alJEoOSi : ISalim et ^l2007h for both red sequence 
and blue cloud galaxies (scatter ~0.2 dex). 

In addition to the stellar masses the absolute magnitudes and 
K-corrections were derived from the output model stellar popula¬ 
tions of this fitting procedure. K-corrections for each galaxy were 
obtained by interpolating the model u,g,r,i,z magnitudes. For this 
work we require the kng = ku — kg correction to obtain the u-g rest 
frame colour and also kr to construct volume limited samples. We 
note that K-corrections are very small for the relatively low red- 
shift sample studied here (0.02 < z < 0.085) and the analysis 
presented in this paper does not critically depend on the accuracy 
in K-correction calculations. 

3.2 Photometric redshifts 

The SDSS database provi des two different types of pho tomet- 
ric redshift measurements (lAdelman-McCarthv et al.|[200^ . One 
of the measurements is based on a template fitting approach 
dCsabai & Budavarill2003h and the other on an artificial neural net- 
work approach dCollister & Lahavl2013h . The fitting approach uses 
a set of galaxy spectral templates with a range of redshifts, types 
and luminosities. Expected colours are derived from the templates 
and fit to the observed galaxy colours. The templates with the best 
fits are found, from which photometric redshifts are obtained. 

The artificial neural network consists of a number of layers 
of nodes with connections between nodes in adjacent layers. The 
network accepts the galaxy photometry and outputs photometric 
redshifts. A representative training set of galaxies is required to 
tune the connections between the nodes. Once the connections have 
been tuned the network can quickly compute photometric redshifts 
from galaxy photometry. 

In this work we show results based on photometric redshifts 
obtained from the template fitting method. Although not presented 

^ http://www.mpa-garching.mpg.de/SDSS/DR7/ 


in this paper we have also checked some results using photomet¬ 
ric redshifts obtained from an artificial neural network method. 
Quantitatively similar results were obtained and hence the choice 
of method does not affect the conclusions. 

For the SDSS photometric redshift sample the derived r-band 
absolute magnitudes based on the spectroscopic redshifts were ad¬ 
justed to account for the difference in distance modulus between 
the spectroscopic and photometric redshift measurements. 

We calculated the redshift uncertainty of the photometric red¬ 
shifts by taking half the difference between the 16^^ and 84^^ per¬ 
centiles of the distribution of the difference between the spectro¬ 
scopic and photometric redshifts. We find a redshift uncertainty of 
0.0185 for the SDSS photometric redshift sample. 

3.3 Simulated photometric redshifts 

In addition to the photometric redshifts from the SDSS photome¬ 
try simulated photometric catalogues were created by displacing 
galaxies from their spectroscopic positions. Displacements were 
obtained for all of the galaxies that had spectroscopic redshift mea¬ 
surements. The redshift displacements were generated by drawing 
them from a Gaussian distribution with a particular standard devi¬ 
ation. The standard deviation (cr^) is the redshift uncertainty of the 
catalogue. Simulated photometric redshifts were created by adding 
the redshift displacements to the spectroscopic redshifts. For the 
simulated photometric redshift samples the derived r-band abso¬ 
lute magnitudes based on spectroscopic redshifts were adjusted to 
account of the difference in distance modulus between the spec¬ 
troscopic and simulated photometric redshift measurements. In to¬ 
tal 14 simulated catalogues were created with the following red¬ 
shift uncertainties: 0.0025, 0.005, 0.0075, 0.01, 0.0125, 0.0150, 
0.0175, 0.02, 0.03, 0.04, 0.05, 0.06, 0.08 and 0.1. 

3.4 Volume limited samples 

Intrinsically faint galaxies become undetectable as redshift in¬ 
creases. An apparent magnitude limited sample would, therefore, 
suffer from a volume bias. To mitigate this effect volume lim¬ 
ited samples were constructed. This was achieved by plotting the 
K-corrected r-band absolute magnitude, Mr, as a function of the 
spectroscopic redshift. Figure [2] shows the selection of the spectro¬ 
scopic redshift sample. The vertical magenta lines mark the lower 
(z = 0.02) and upper (z = 0.085) bounds of the sample. We com¬ 
pute the r-band absolute magnitude detection limit using the appar¬ 
ent magnitude limit of the survey (17.77), the distance modulus and 
the 95* percentile of the r-band K-corrections. The magenta con¬ 
necting curve in Figure |2] shows this detection limit. The limiting 
magnitude, Mr (limit), at the upper redshift boundary (z = 0.085) 
was found to be —20.26. Selecting only those galaxies brighter 
than this value gives the volume limited sample. The same limit¬ 
ing magnitude obtained for the spectroscopic catalogue was then 
applied to the photometric and simulated catalogues. The galaxies 
in the volume limited samples are used to trace the density field and 
are hence also the density defining population. 

3.5 Target sampling rate 

The observational setup used in the SDSS consists of optical fibers 
mounted on aluminium plates. The finite size of the fibers makes 
it difficult to target more than one galaxy within an angular radius 
of 55 arcsecs using a single plate. Targeting two galaxies separated 
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Table 1. Environment methods and parameters 



Method 

Parameters 

(i) 

A/^* nearest neighbour 

A=l-10; du=1000-20000 km/s 

(ii) 

Fixed aperture 

r=0.1-3.0 Mpc; du=1000-20000 km/s 


by less than this would result in a fiber collision. This targeting 
results in systematic under sampling that depends on environment. 

A target sampling correction is required to account for t his. _ 

_ ^Following previous work (iBaldrv & Balogh|[2006l : IPeng et aP 

l201Qh . the target sampling rate was calculated for each target in the 
spectroscopic redshift, photometric redshift and simulated photo¬ 
metric redshift volume limited samples. The target sampling rate 
is obtained by dividing the number of galaxies within the volume 
limited sample and within 55 arcsec of the target by the number 
of galaxies in the photometric sample, including galaxies with no 
redshift but would have bright enough absolute magnitudes to enter 
the volume limited sample within our redshift range, and within 55 
arcsec of the target. A target sampling rate correction is required 
for ~ 10 percent of the objects. The target sampling rate correction 
was used in two places: (i) to weight (x l/(target sampling rate)) the 
contribution of each galaxy in the computation of galaxy environ¬ 
ment. (ii) to weight (x l/(target sampling rate)) each galaxy when 
representing the overall galaxy population. 


3.6 Environment methods 

In this work we calculate a number density (pi) for each galaxy 
in the photometric, spectroscopic and simulated samples using two 
methods: (i) the nearest neighbour method and (ii) the fixed 
aperture method. Both methods have two aperture parameters: one 
to constrain the projected size of the aperture on the sky and the 
other to constrain the depth of the aperture (see TableB. 

In the nearest neighbour method velocity cuts (±dv km/s) 
centered on the target galaxy are used to constrain possible neigh¬ 
bours. The Hubble fiow velocity is assumed to be the speed of light 
multiplied by the redshift. The angle between the target and these 
galaxies are then computed. The galaxy with the smallest an¬ 
gle is identified. This is the A^* nearest neighbour. This angle is 
used to define the aperture volume which is approximately a coni¬ 
cal frustum centred on the target. The number density is then sim¬ 
ply the sum of the N target sampling corrections divided by the 
volume of the aperture. In this study we adopt N = 1 — 10 and 
dv = 1000 — 20, 000 km/s. For galaxies close to the redshift 
boundaries we adjust the aperture volumes as described in Section 

Ezl 

In the fixed aperture method again velocity cuts (±dv km/s) 
centered on the target galaxy are used to constrain the depth of the 
aperture. The projected size is fixed by choosing a value for the 
radius (r Mpc) of the aperture at the redshift of the target galaxy. 
Again the aperture is approximately a conical frustum. We use r = 
0.1 —3.0 Mpc and du = 1000 — 20, 000 km/s. For galaxies close to 
the redshift boundaries we adjust the aperture volumes as described 
in Section im Fixed apertures that are devoid of density defining 
galaxies are assigned a nominal minimum density of 0.5 galaxies 
per aperture. 


After we have computed a density for each galaxy we then 
evaluated an overdensity, (5, using the equation below: 

d — (^Pi Pm) jPm (1) 

The number density, pi, is the density described above calcu¬ 
lated using either the A/* nearest neighbour or fixed aperture meth¬ 
ods. Pm is the mean density of galaxies (weighted by the target 
sampling rate corrections) within a redshift window centered on 
the target galaxy utilizing all of the available area. Typically the 
size of the redshift window was set to be the same scale as the 
velocity cuts used by the method. We finally compute the galaxy 
environment using the expression: 

Lo^l0(l + (5) (2) 


3.7 Survey Edge 

The SDSS footprint is not a simple geometry in right ascension 
and declination. It consists of a series of parabolic strips. The 
edge of the footprint is irregular (see Figure [TJ- To manage the 
edges of the survey often a survey mask populated with a set of 
random galaxies is employed. This is especially common in sce- 
narios where the two-point correlation function is required (e.g. 
IWang et ^l2013h . For simplicity, however, we constructed a high 
order polygon (> 1000) within the footprint to define the edge. 
This polygon has an area of 6926.64 square degrees. Galaxies out¬ 
side of this perimeter are discarded. In total ~10 percent of the 
galaxies were discarded but many of these galaxies were in stripes 
outside of the main footprint. The angular separations between 
galaxies inside the region and the closest points on the perimeter 
of the polygon were computed. 

For a galaxy close to the perimeter a fraction of the aper¬ 
ture volume used to calculate the environment would reside out¬ 
side of the survey where there is no data. Environments computed 
for galaxies at the edges of the survey footprint would, therefore, 
tend to be underestimated. To manage this, the fraction, /, defined 
as the distance to the edge of the polygon divided by the projected 
radius of the aperture was calculated. We adopt a conservative ap¬ 
proach and discard targets with / < 1 . 0. Oth er authors use simi¬ 
lar approaches. For example IPeng et al.l fcOlOh consider only those 
galaxies where / > 0.9 and applies a correction to the aperture vol¬ 
umes where / > 0.9 but / < 1.0. Another approach (ICooper et al.l 
l2005h is to discard galaxies that are within a certain distance from 
the perimeter (e.g. 1 Mpc). 

In addition to the angular boundaries we manage the redshift 
boundaries. In cases where the target galaxies are close to a redshift 
boundary (i.e. the aperture would cross over the boundary) we re¬ 
duce the depth of the half of the aperture infringing the boundary. 
We fix the half depth of this half of the aperture to be the comoving 
distance from the target galaxy to the boundary. This ensures the 
aperture fits inside the redshift range. The depth of the other half of 
the aperture that resides within the redshift range is not altered. 


4 RESULTS 

We present the results in four sections. In Section l4Jl we compare 
the scales probed by the A/* nearest neighbour and fixed aperture 
methods and establish spectroscopic benchmark environments. In 
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Figure 3. Shows the median comoving distance to the nearest neigh¬ 
bour for a range of values of AT as a function of the velocity cut {dv). 
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Log(l+(5) V*'’ Nearest Neighbour 

Figure 4. Shows the hxed aperture (r = 1.8 Mpc and dv = 1000 km/s) 
spectroscopic environment measurements versus the nearest neighbour 
(AT = 4 and dv = 1000 km/s) spectroscopic environment measurements. 
The number of galaxies in each 2-dimensional bin is represented with a 
colour as indicated in the colour bar. The linear best ht is shown with the 
solid black line and the 1-to-l is shown with the dashed line. The correlation 
quoted is Spearman rank correlation coefficient. 

Section 1431 we compare the spectroscopic environment measure¬ 
ments from A/^* nearest neighbour and fixed aperture methods. In 
Section 14.31 we compare the spectroscopic and photometric envi¬ 
ronment measurements for the two methods. In Section |4!41 we in¬ 
vestigate how the Spearman rank correlation coefficient between 
environment measurements and the spectroscopic benchmark en¬ 
vironment measurements varies with the aperture parameters and 
redshift uncertainty. We determine the aperture parameters that give 
the strongest correlations as a function of redshift uncertainty. 


4.1 Spectroscopic benchmark environments 

To establish spectroscopic benchmark environments to facilitate 
comparisons between environment measurements as a function of 


the aperture parameter values and redshift uncertainty we examine 
the comoving scales that the methods probe. 

In both methods the depth of the aperture is controlled by 
the same parameter: the velocity cut. In this study we consider a 
range of velocity cuts: 1000 — 20, 000 km/s. A velocity cut of 
20, 000 km/s gives an aperture that extends through the entire red¬ 
shift range: 0.02 — 0.085. The possible comoving aperture depths 
at the median redshift are: 27.9 — 271.7 Mpc. 

Next we examine the size of the aperture perpendicular to the 
line of sight for the A/^* nearest neighbour method. Figure [3 shows 
the median A^* nearest neighbour distance as a function of the ve¬ 
locity cut using spectroscopic measurements for values of N in the 
range: 1 — 10. The plot illustrates the obvious points that as N 
increases the median distance to the A^* neighbour also increases 
and as the velocity cut increases the median distance to the A7* 
nearest neighbour decreases. More importantly Figure [3 shows that 
the median A/^* nearest neighbour distance for this range of N is 
~0.3 — 3.4 Mpc. In this study we choose to calculate the A7* near¬ 
est neighbour environments using N — 1 — 10 and fixed aperture 
environments using r = 0.1 — 3.0 Mpc so that the range of scales 
probed by the two methods are matched. 

We establish two sets of spectroscopic benchmark environ¬ 
ments: one set for the fixed aperture method and one set for the 
A^* nearest neighbour method. To strike a balance between probing 
small scales ( ^1.0 Mpc) where environment effects ye rep orted to 
be important (iBlanton & Berlindll2007l : IWilman et al.l201Ql) . ensur¬ 
ing we have a large dynamic range of environments and we do not 
have too many apertures devoid of galaxies we choose an aperture 
scale of 1.8 Mpc perpendicular to the line of sight and a velocity 
cut of 1000 km/s for the fixed aperture spectroscopic benchmark 
measurements. Referring to Figure [3 to ensure that the aperture 
scales for the two methods are well matched we choose the aper¬ 
ture parameters A^ = 4 and dv = 1000 km/s for the A^* nearest 
neighbour spectroscopic benchmark measurements. 

4.2 Comparison of spectroscopic environments 

Figure|4]shows the A^* nearest neighbour spectroscopic benchmark 
environments versus the fixed aperture benchmark spectroscopic 
environments. N = 4 and dv = 1000 km/s are adopted for the A/^* 
nearest neighbour method and r = 1.8 Mpc and dv = 1000 km/s 
are adopted for the fixed aperture method. The dynamic range for 
the fixed aperture method is smaller than the A^* nearest neigh¬ 
bour method and so the best linear fit (gradient of 0.56) relation¬ 
ship shown with the solid black line deviates from the 1-to-l line 
shown with the dashed black symbols. Although we show a linear 
fit on the plot the relationship between the environments is non¬ 
linear. For large and small environments the A^* nearest neighbour 
method tends to ‘stretch out’ the fixed aperture environments to 
higher and lower values respectively. Nevertheless there is a strong 
correlation between the spectroscopic environments measured with 
the two methods. The Spearman rank correlation coefficient is 0.91. 

4.3 Spectroscopic vs photometric environments 

Figure |5] shows density plots of the spectroscopic benchmark envi¬ 
ronments versus the SDSS photometric environments on the left 
and the spectroscopic (red) and photometric (blue) environment 
distributions on the right for the A^* nearest neighbour (top) and 
fixed aperture method (bottom). The benchmark parameters are 
adopted (A^ = 4 and dv = 1000 km/s for the A^* nearest neigh¬ 
bour method and r = 1.8 Mpc and dv = 1000 km/s for the fixed 
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Figure 5. Density plots of the spectroscopic environment measurements versus the photometric environment measurements for the N'-^ nearest neighbour 
method (a) and the hxed aperture method (c). The corresponding environment distributions for the spectroscopic (red) and photometric (blue) measurements 
are shown in plots (b) and (d) for the AA* nearest neighbour method and the fixed aperture method respectively. The aperture parameters used for the 
spectroscopic measurements are N = 4 and dv = 1000 km/s for the nearest neighbour method and r = 1.8 Mpc and dv = 1000 km/s for the fixed 
aperture method. The aperture parameters for the photometric measurements are A/^ = 7 and dv = 7000 km/s for the A^* nearest neighbour method and 
r = 1.9 Mpc and dv = 6000 km/s for the fixed aperture method. The linear best fit lines for the measurements are shown in black and 1-to-l lines are dashed. 
The correlations quoted are Spearman rank correlation coefficients. 


aperture method) for the spectroseopie measurements whereas we 
choose N = 7 and dv = 7000 km/s for the A7* nearest neigh¬ 
bour method and r = 1.9 Mpc and dv = 6000 km/s for the fixed 
aperture method for the photometric measurements as we find these 
parameters lead to the strongest correlations (see Section [Jd]). In¬ 
cluded on the density plots are the 1-to-l lines (dashed) and the 
linear best fit lines (solid). The colour of each grid cell on the den¬ 
sity plots represents the number of galaxies in that cell. 

Plots a) and b) show the relatively smooth and continuous 
nature of the A7* nearest neighbour environments. Conversely 
the underdense, fixed aperture, spectroscopic environments shown 
in plots c) and d) are noticeably discretized. The fixed aperture 
method produces forbidden environments because this method is 
based on counting the number of whole galaxies within a volume. 
The spectroscopic fixed aperture distribution does not fall to zero in 


these forbidden ranges, however, for two reasons: (i) ~10 percent 
of the galaxies are weighted by the target sampling rate correction 
and (ii) densities are measured relative to the mean densities at the 
redshifts of the target galaxies. This discretization effect is unno- 
ticeable in the photometric fixed aperture environment distribution 
because a larger velocity cut was adopted leading to more galaxies 
on average being found within the apertures. 

Plots b) and d) show that the photometric environment dis¬ 
tributions are more peaked and their ranges are smaller than the 
corresponding spectroscopic distributions. Quantitatively the stan¬ 
dard deviations of the environment distributions also shows this 
is the case. The standard deviation of the A^* nearest neighbour 
method for the spectroscopic benchmark environment distribution 
is 0.7 whereas the standard deviation for the photometric environ¬ 
ment distribution is 0.4. The standard deviation for the fixed aper- 
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Figure 6. Environment of two target galaxies: one in a low density environment (top row) and one in a high density environment (bottom row) measured with 
spectroscopy (left column) and photometry (right column). The targets are marked in red. The blue circles mark a comoving distance of 1.8 Mpc from the 
target. The blue symbols mark galaxies that are inside the circle. The black symbols mark galaxies that are outside the circle. Only those galaxies that are in 
the density defining population and are also constrained by a velocity cut of 1000 km/s are show. The number of galaxies found within 1.8 Mpc of the target 
for plots a) to d) are: 1, 5, 24 and 13. 


a) 


b) 



Figure 7. Shows the redshift and declination of galaxies within a wedge of volume using spectroscopic measurements (a) and photometric measurements (b). 
The redshift range is: 0.02 — 0.085 and extends along the radial direction of the cones. The declination range is 20 — 40 degrees. The right ascension is 
restricted to the range 215 — 220 degrees. 


ture method reduces from 0.44 to 0.31. The most dense environ¬ 
ments measured with spectroscopy are found to be less dense with 
photometry. This is because imprecise redshifts tend to have the 
effect of moving galaxies away from the centres of dense regions. 
Conversely, low density regions in the spectroscopic sample tend to 
be contaminated, in the photometric sample, by galaxies that have 
scattered into them from higher or lower redshifts. 

We illustrate this with Figures [6| and [T] Figure [6] shows two 
typical target (marked in red) galaxies: one in a low density region 
(top row) and one in a high density region (bottom row). The red¬ 
shifts of the galaxies are measured with spectroscopy in the plots 


on the left and photometry in the plots on the right. The plots show 
the galaxies that are constrained by velocity cuts of ±1000 km/s. 
The low density environment that was measured with spectroscopy 
is measured to be more dense with photometry. The number of den¬ 
sity defining galaxies found within 1.8 Mpc of the target (marked 
by the blue circle) increases from 1 to 5. The high density envi¬ 
ronment that was measured with spectroscopy is measured to be 
less dense when measured with photometry. The number of den¬ 
sity defining galaxies found within 1.8 Mpc of the target decreases 
from 24 to 13. 

Figure |7] shows the redshift and declination of galaxies within 
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Figure 8. Shows the Spearman rank correlation coefficients between the nearest neighbour environments and the spectroscopic benchmark environments 
(where N = A and dv = 1000 km/s) for nine catalogues of galaxies as a function of the aperture parameters: N and the velocity cut (dv). a) Shows the 
correlations for the spectroscopic redshift catalogue, b) Shows the correlations for the photometric redshift catalogue, c) to i) show the correlations for simulated 
photometric redshift catalogues with redshift uncertainties of 0.0025, 0.005, 0.01, 0.02, 0.04, 0.06 and 0.1 respectively. The correlation is represented with a 
colour in each grid cell. The equally spaced contours in the logarithm of the correlation highlight the location of the peak in the parameter space on each grid. 


a wedge of volume using spectroscopic measurements (a) and pho¬ 
tometric measurements (b). Clustering can easily be seen in the 
spectroscopic plot whereas in the photometric plot the distribution 
of galaxies looks more random. It is difficult to see the clustering 
(from this perspective) in the photometric plot because of the red¬ 
shift errors. High density regions are smoothed into the surrounding 
low density regions leading to a more homogenous density field. 

In Figure[5]we also computed the gradients of the linear best fit 
lines and Spearman rank correlation coefficients for plots a) and c). 


The gradient of the linear best fit lines between the spectroscopic 
and photometric environment measurements are 0.38 for the 
nearest neighbour method and 0.45 for the fixed aperture method. 
We note that although we have shown linear best fit lines on the 
plots the relationships between the spectroscopic and photometric 
environments are nonlinear. The high density photometric environ¬ 
ments for both methods deviate away from the linear fits to larger 
values. We find in agreement with previous studies that the high 
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Figure 9. Shows the Spearman rank correlation coefficients between the hxed aperture environments and the spectroscopic benchmark environments (where 
r = l.S and dv = 1000 km/s) for nine catalogues of galaxies as a function of the aperture parameters: r and the velocity cut (dv). a) Shows the correlations for 
the spectroscopic redshift catalogue, b) Shows the correlations for the photometric redshift catalogue, c) to i) show the correlations for simulated photometric 
redshift catalogues with redshift uncertainties of 0.0025, 0.005, 0.01, 0.02, 0.04, 0.06 and 0.1 respectively. The correlation is represented with a colour in 
each grid cell. The equally spaced contours in the logarithm of the correlation highlight the location of the peak in the parameter space on each grid. 


density environ ment are recovered better than the low density en¬ 
vironments (e.g lCooper et al.ll2005h . 

The Spearman rank correlation coefficients for the nearest 
neighbour method are 0.60 and 0.62 for the fixed aperture method. 
These large correlations suggest that, although it might not be pos¬ 
sible to use photometry to measure galaxy environment for individ¬ 
ual galaxies, with a high degree of accuracy, it should be possible 
to extract an environment signal from a sufficiently large sample of 
galaxies. 


4.4 Correlation vs aperture parameters and redshift 
uncertainty 

Figure [8] and Figure |9] show how the Spearman rank correlation co¬ 
efficient between environments and the spectroscopic benchmark 
environments vary as a function of the aperture parameters and red¬ 
shift uncertainty for the nearest neighbour and fixed aperture 
methods respectively. Each figure consists of a sequence of nine 
coloured grids. Each grid shows the results for a different redshift 
uncertainty. The results for the SDSS spectroscopic redshift sample 
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Figure 10. Left column shows the optimal values for the parameters that control the projected size of the apertures as a function of redshift uncertainty. 
Middle column shows the optimal velocity cut as a function of redshift uncertainty. Right column shows the Spearman rank correlation coefficient between 
the environments and the spectroscopic benchmark environments achieved using the optimal parameter values. The top row is for the nearest neighbour 
method and the bottom row is for the fixed aperture method. The results for the SDSS photometric redshift are shown with the star symbols. 


are shown in plot a). The results for the SDSS photometric redshift 
sample are shown in plot b) and the results for the simulated pho¬ 
tometric samples with az — 0.0025, 0.005, 0.01, 0.02, 0.04, 0.06 
and 0.1 are shown in plots c) to i) respectively. In Figure[8]the aper¬ 
ture parameter: N varies along the vertical axis whereas in Figure 
mthe radius of the fixed aperture varies along the vertical axis. The 
velocity cut (dv) used to constrain the depth of the apertures varies 
along the horizontal axis of each grid. The colour of each grid cell 
represents the Spearman rank correlation coefficient between the 
environments and the spectroscopic benchmark environments. The 
smallest correlations are shown in blue and the largest correlations 
are shown in red. The equally spaced contours in the logarithm of 
the correlation highlight the location of the peak in the parameter 
space on each grid. 

From visual inspection of Figure [8] and Figure [9] the correla¬ 
tion between environments and the spectroscopic benchmark en¬ 
vironments for the fixed aperture method and nearest neigh¬ 
bour method behave in a similar way. As the redshift uncertainty 
increases the coloured grids transition from red to blue. The Spear¬ 
man rank correlation coefficient between the environments and the 
spectroscopic benchmark environments therefore decreases as the 
redshift uncertainty increases. 

The Spearman rank correlation coefficient between the envi¬ 
ronments and the spectroscopic benchmark environments depends 
not only on the redshift uncertainty but also on the aperture pa¬ 
rameters. For low redshift uncertainties (az < 0.02) there exists a 
‘sweet spot’ (e.g. see plot e) in either Figure). A prudent choice 
for the aperture parameter values enables the strongest environ¬ 
ment signal to be extracted. At larger redshift uncertainties pro¬ 
viding the aperture parameter values are large enough (e.g. N > 3, 
r > 1.2 Mpc, dv > 10, 000 km/s and smaller than the maximum 


values we study) the correlation is relatively insensitive to the aper¬ 
ture parameter values. 

Figure [T0| shows the parameter values that give the largest 
Spearman rank correlation coefficient between the environments 
and spectroscopic benchmark environments as a function of red¬ 
shift uncertainty. The figure consists of six plots. The top row (a, b 
and c) are for the A^* nearest neighbour method and bottom row (d, 
e and f) are for the fixed aperture method. The left hand column (a 
and d) shows the optimal parameter values for the parameters that 
control the projected sizes of the apertures. The best values of N 
are shown in plot a) and the best values of r are shown in plot d) 
as a function of redshift uncertainty. The middle column (b and e) 
shows the optimal velocity cut as a function of redshift uncertainty 
for the two methods and the right hand column (c and e) shows 
the the Spearman rank correlation coefficient obtained with the op¬ 
timal parameter values as a function of redshift uncertainty. The 
results for the SDSS photometric redshifts are shown with the star 
symbols. 

The optimal parameter values for the parameters that control 
the projected size of the apertures are aligned with the values cho¬ 
sen for the spectroscopic benchmark measurements. Both N and 
r are fairly insensitive to the redshift uncertainty (there appears to 
be a mild trend with increasing redshift uncertainty). N = 5 and 
r = 1.8 Mpc are good choices across the entire range of redshift 
uncertainties. 

The optimal velocity cut for both methods however is depen¬ 
dent on the redshift uncertainty. To obtain the optimal correlation 
with the spectroscopic benchmark measurements it appears that the 
A^* nearest neighbour method requires a slightly deeper aperture 
than the fixed aperture method. Nevertheless for both methods as 
the redshift uncertainty increases a larger velocity cut should be 
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employed. This ensures that galaxies displaced from their spectro¬ 
scopic positions along the line of sight by uncertain redshift mea¬ 
surements are captured by the apertures. 

The optimal parameter values and correlations obtained using 
the SDSS photometric redshift sample (which has an uncertainty of 
0.0185) are consistent with the values obtained with the simulated 
photometric redshifts. 

The best Spearman rank correlation coefficient between the 
environments and the spectroscopic benchmark environments de¬ 
creases from 1.0 (i.e. when comparing the benchmark environ¬ 
ments with themselves) down to ~0.4 with simulated redshifts with 
an uncertainty of 0.1 for both the nearest neighbour and fixed 
aperture methods. The correlation decays rapidly for small red¬ 
shift uncertainties {gz < 0.02) and then more gradually for larger 
redshift uncertainties {gz > 0.02). Although the correlation at a 
redshift uncertainty of 0.1 is reduced it is significantly larger than 
zero. It therefore should be possible to extract an environment sig¬ 
nal from large scale photometric surveys such as DES. 


5 DISCUSSION 

This section consists of three parts. In Section [5Jl we examine the 
environmental dependence of the galaxy red fraction using pho¬ 
tometric redshifts. In Section 15.21 we estimate how much larger, 
than spectroscopic samples, photometric samples must be to ob¬ 
tain equivalent measurements of environment correlations and in 
Section [53] we briefiy compare this work with others in the litera¬ 
ture. 

5.1 Red fraction: mass-environment dependence 

In the local universe it has been shown that the fra ction of red galax¬ 
ies d e pends on stellar mass and environment (iBaldrv & BaloghI 
bood IPengeTani2010h . As a case study we investigate the de¬ 
pendence of the galaxy red fraction on environment measured with 
phot ometric red s hifts. 

jPeng et akl fcOlOh showed that in the local universe the red 
fraction is high for massive galaxies even in low density envi¬ 
ronments and the red fraction is high for galaxies in very dense 
environments even for galaxies with low masses. Recent work 
has shown that th is dependence is establish at least at z = 0.7 
(iKovac et al .1120141) . 

These results are consistent with the idea of two main chan¬ 
nels to quench star formation: “mass quenching” and “environment 
quenching”. These quenching modes act independently and facili¬ 
tate the transition of galaxies from the “blue cloud” to the “red se¬ 
quence”. One candidate for “mass quenching” is AGN feedback. 
In this scenario the central black hole drives energetic outfiows 
which can heat gas and hence inhibit star formation. “Environment 
quenching” could manifest itself when satellite galaxies fall into 
larger dark matter haloes. Ram pressure stripping or strangulation 
are physical mechanisms that could be responsible for this quench¬ 
ing mode. 

To study the dependence of red fraction on mass and environ¬ 
ment as a function of redshift uncertainty, here, we extend the range 
of masses down to 10®'^ M©. To do this we calculate environments 
for galaxies that are intrinsically fainter than the absolute magni¬ 
tude cut off used to construct the volume limited samples previ¬ 
ously. These faint galaxies however do not become part of the den¬ 
sity defining populations. To account for the fact that these galaxies 
are not detectable through the entire redshift range we determine a 


volume correction for these galaxy. The volume correction is the 
ratio of the total volume imposed by the redshift range to the de¬ 
tectable volume of the galaxy with a particular intrinsic brightness. 
The volume correction is used to weight galaxies (together with 
the target sampling correction) when representing the whole popu¬ 
lation of galaxies. 

We now outline how we determined the red fraction for the 
galaxies in the samples. Eor each galaxy the SDSS apparent magni¬ 
tudes: u,g,r,i and z were adjusted to the AB system using small cal - 
ibration offsets (lEisenstein et al.|[200^ : [Blanton & Row^l2007l) . 
The (u-g)Rest colour for each galaxy was calculated by subtracting 
the K-correction (kug) obtained from the stellar population fits. 

(u - ^)Rest = (u - K) - {g- kg) = {u- g) - kug (3) 

The (U-B)Rest colour was then calculated using the transfor¬ 
mation equation by Lupton (2005) found on the SDSS website: 

{U - B)Rest = 0.8116((u - ^)Rest) - 0.1313 (4) 

The colour-mass diagrams for the spectroscopic and photo¬ 
metric selections are shown in Eigure [TT] The bimodality in the 
galaxy population can clearly be seen in both samples. Eollowing 
[Peng et aP feOlOh a dividing line was employed, restated here: 

{U - B)Rest = 1.10 + 0.075/o^(m/10^°Mo) - 0.18^ (5) 

The galaxies with a (U-B)Rest colour larger than the dividing 
value, for the mass and redshift, defined by equation 0 were con¬ 
sidered to be red and the galaxies with a (U-B)Rest colour smaller 
than the dividing value were marked as blue. This binary division 
facilitated a simple way to compute the fraction of red galaxies in 
any bin of mass and environment. The red fraction we used was 
simply the number of red galaxies, weighted by the target sampling 
rate and volume corrections, in the bin divided by the total weighted 
number of galaxies in the bin. 

Eigure [12] shows plots of the galaxy red fraction as a function 
of log(mass) and environment for the fixed aperture method (left 
column) and the nearest neighbour method (middle and right 
columns). In the left and middle columns the derived masses are 
based on the spectroscopic redshifts. In the right column the derived 
masses are also based on the spectroscopic redshifts but the masses 
have been convolved with a Gaussian distribution with a width of 
0.3 dex. This is to account for e rrors in the mass estimates due 
to the photometric redshifts fe.g. [Tavlor et alJ[2009[ : [Eossati et alJ 
[2015[) . The SDSS spectroscopic redshift sample is shown in the 
first row (a and b), the SDSS photometric redshift sample is shown 
in the second row (c, d and e) and simulated photometric redshift 
samples with redshift uncertainties of 0.4 and 0.6 are shown in the 
third (f, g, h) and fourth (i, j, k) rows respectively. The aperture 
parameter values used are close to the optimal values determined 
in Section [43] and are stated on each plot. The contours mark the 
red fractions: 0.2, 0.4, 0.6 and 0.8 

The first row of plots (a and b) for the spectroscopic redshift 
sample show that the the red fraction increases with both log(mass) 
and environment. There are clear contours of constant red fraction 
covering several orders of magnitude in log(mass) and environment 
for both the A* nearest neighbour method (b) and the fixed aper¬ 
ture m ethod (a). We reproduce the result presented bv [Peng et al.[ 
boioh using spectroscopic redshifts using the A* nearest neigh¬ 
bour method and in addition the fixed aperture method. The red 
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Figure 11. (U-B) rest frame colour as a function of the log(mass) for the spectroscopic selection (left) and the photometric selection (right). The colour 
represents the relative density of the galaxies in the parameter space with the target sampling rate and volume corrections applied. Beneath a threshold the 
colour scheme is not used but the individual galaxies are overplotted with the “+” symbol. The black dividing line (specihed in the main text) splits the red 
sequence and blue cloud populations. 


fraction is high for massive galaxies even in low density environ¬ 
ments and the effect of environment is most important for galaxies 
with low masses. 

The second row of plots (c, d and e) for the SDSS photometric 
redshift sample also show contours of constant red fraction, albeit 
some deterioration. The information to separate the red fraction by 
log(mass) and environment is still present in the SDSS photomet¬ 
ric redshift measurements. Galaxies with smaller masses than in 
the spectroscopic sample are found in the photometric samples be¬ 
cause some galaxies with small spectroscopic redshifts (z < 0.02) 
are scattered into the sample because of the uncertain redshift mea¬ 
surements. 

As noted earlier (in Section 14.31 1 for fixed aperture parame¬ 
ters the range of environments for the photometric measurements 
are smaller than the range for the spectroscopic measurements. The 
reduction in dynamic range is particularly evident when compar¬ 
ing the red fraction surfaces for the SDSS spectroscopic and pho¬ 
tometric redshift samples for the nearest neighbour method 
(plots b and d). The red fraction surface for the SDSS photomet¬ 
ric redshift sample is shrunk along the environment axis compared 
with the spectroscopic redshift sample. The dynamic ranges for the 
SDSS photometric redshift and spectroscopic redshift samples for 
the fixed aperture method (plots a and c) are however similar. In 
this method the much larger velocity cut (x6) adopted for the SDSS 
photometric redshift (plot c) is able to compensate for the reduction 
in dynamic range due to the photometric redshifts. 

Increasing the redshift uncertainty to 0.4 and 0.6 (third and 
fourth rows) results in more deteriorated red fraction surfaces and 
smaller dynamic ranges. Nevertheless, particularly for the 
nearest neighbour method the red fraction contours still behave in a 
similar manner to the SDSS spectroscopic sample. The red fraction 
dependence on galaxy environment does not break down even for 
samples with large redshift uncertainties. The A^* nearest neigh¬ 
bour method appears to fare better than the fixed aperture method 
as the redshift uncertainty increases. This is because the A^* nearest 
neighbour method tends to probe a larger dynamic range of envi¬ 
ronments than the fixed aperture method. 

To model the impact of the photometric redshifts on the mass 
estimates we have applied an additional uncertainty of 0.3 dex to 


the mass estimates in the right hand column. Both the middle and 
right hand columns employ the A/^* nearest neighbour method so 
comparing the red fraction surfaces in these columns shows the im¬ 
pact of the additional mass uncertainty. The additional uncertainty 
in the mass estimates tends to mix up the red and blue galaxies re¬ 
sulting in more homogenous red fractions. However even with this 
additional uncertainty the general trends we described above are 
still present. 


5.2 Photometric sample size for environment correlations 

One of the main goals of this paper is to quantify the impact of 
uncertain redshifts on measurements of environment trends. The 
questions we consider here are: (i) How many times larger than the 
spectroscopic sample must the photometric sample be to measure 
environment correlations to the same accuracy? (ii) How does this 
factor depend on the trend correlation with spectroscopic measure¬ 
ments? (iii) How does this factor vary with redshift uncertainty? 

To calibrate our measurements to enable comparisons between 
the results obtained with different methods, aperture parameter 
values and redshift uncertainties we constructed two sets of fake 
data which we label: F-Low and F-High. The fake data represents 
some galaxy property that increases monotonically with environ¬ 
ment. More specifically the fake properties were constructed to be 
linearly related to both the A^* nearest neighbour spectroscopic 
benchmark environments (N = 4 and dv = 1000 km/s) and 
the fixed aperture spectroscopic benchmark environments (r = 
1.8 Mpc and dv = 1000 km/s) with a fixed level of Gaussian 
scatter. F-Low was constructed to have a Spearman rank correla¬ 
tion coefficient of 0.3 whereas F-High was constructed to have a 
Spearman rank correlation coefficient of 0.7 with both sets of spec¬ 
troscopic benchmark environments. Two important simplification 
to note are that the value of the fake property does not change with 
the value of galaxy’s redshift and no mass dependence is built into 
the fake properties. Figure[T3shows plots of the fake properties as a 
function of the spectroscopic benchmark environments for the A/* 
nearest neighbour (a and c) and fixed aperture (b and d) methods. 

We calculated the Spearman rank correlation coefficient be¬ 
tween the environments and the fake properties as a function of 
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Figure 12. Shows the red fraction as a function of Log(mass) and environment for the fixed aperture method (left column) and the nearest neighbour 
method (middle and right columns). The derived masses are based on the spectroscopic redshifts in the left and middle columns. The masses of the galaxies in 
the plots in the right hand column have been convolved with a Gaussian distribution with a width of 0.3 dex to model the impact of the photometric redshifts 
on the mass estimates. The spectroscopic redshift sample is shown in the first row (a and b), the photometric redshift sample is shown in the second row (c, d 
and e) and simulated photometric redshift samples with redshift uncertainties of 0.04 and 0.06 are shown in the third (f, g and h) and fourth rows (i, j and k). 
The aperture parameter values used are stated on each plot. The contours mark the red fractions: 0.2, 0.4, 0.6 and 0.8. 
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Figure 13. Shows the spectroscopic benchmark environment measurements for the nearest neighbour method (a and c) and the fixed aperture method (b 
and d) versus two sets of faked data. One set of faked data was constructed to have a low correlation (a and b) and the other a high correlation (c and d) with 
the spectroscopic benchmark environments. The aperture parameters for the spectroscopic benchmark measurements are A/' = 4 and dv = 1000 km/s for the 
A^* nearest neighbour method and r = l.S Mpc and dv = 1000 km/s for the fixed aperture method. 


redshift uncertainty and the aperture parameters. We made figures 
similar to Figures [8] and |9] and then found the aperture parameters 
that result in the largest correlations for each redshift uncertainty. 
The parameters for each redshift uncertainty are shown in Table |2] 
The second and third columns of the table show the parameters for 
the A^* nearest neighbour method. The fourth and fifth columns of 
the table show the parameters for the fixed aperture method. 

Next we generated 224 different realizations of the redshift 
catalogues. To obtain different realizations of the spectroscopic cat- 
alogue we assigned a nominal redshift uncertainty of 0.0001 (e.g. 
ICooper et al.ll2005h to the spectroscopic redshifts. For each realiza¬ 
tion we followed the same procedure as before but used the param¬ 
eters specified in Table |2] We calculated the Spearman rank cor¬ 
relation coefficient between the environments and the fake galaxy 
properties as a function of the number of objects in the samples. 

Figure [14] shows histograms of the Spearman rank correlation 
coefficients between environment and F-Low for 224 realizations 


of the spectroscopic catalogue (red) and the simulated photometric 
catalogue (blue) with cr^ = 0.02. In this example we calculated 
the Spearman rank correlation coefficient using 10, 000 randomly 
selected objects from each realization. The plot shows that a larger 
redshift uncertainty leads to a smaller Spearman rank correlation 
coefficient. The median correlations are 0.30 and 0.70 for the spec¬ 
troscopic realizations and 0.19 and 0.45 for the photometric real¬ 
izations (with az = 0.02) for F-Low and F-High respectively. In 
addition the correlation distribution is slightly wider for the pho¬ 
tometric realizations. The standard deviation of the Spearman rank 
correlation coefficients of the spectroscopic realizations are 0.0085 
and 0.005 and the standard deviation of the photometric realiza¬ 
tions are 0.0093 and 0.0085 for F-Low and F-High respectively. 
We compute the fractional error for each set of realizations by tak¬ 
ing the ratio of standard deviation of the distribution and the median 
of the distribution. The fractional errors for the spectroscopic real¬ 
izations are 0.029 and 0.0072 and for the photometric realizations 
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Figure 14. Histograms of the Spearman rank correlation coefficients be¬ 
tween the optimal nearest neighbour environments and the faked data 
(constructed to have a correlation of 0.3 with the spectroscopic benchmark 
environments) for 224 realizations of the spectroscopic catalogue (red) with 
(Tz = 0.0001 and 224 realizations of the photometric catalogue (blue) with 
(7z = 0.02. The Spearman rank correlation coefficient is calculated using 
10, 000 objects from each realization. 


Table 2. Aperture parameters to optimize the correlation between the faked data 
and environment for the nearest neighbour method and the fixed aperture 
method for a range of redshift uncertainties. 


Gz 

N 

dv (km/s) 

r (Mpc) 

dv (km/s) 

0.0001 

4 

1000 

1.8 

1000 

0.02 

5 

11000 

1.6 

8000 

0.04 

7 

13000 

1.7 

13000 

0.06 

6 

12000 

2.0 

14000 

0.08 

6 

16000 

2.0 

17000 

0.1 

9 

20000 

2.0 

20000 


The second and third columns are the optimal parameter values for the 

nearest neighbour method and the fourth and fifth colums are the 
optimal parameter values for the fixed aperture method. 

the fractional errors are 0.048 and 0.019 for F-Low and F-High 
respectively. 

Figure [TSl shows four plots of the fractional error as a function 
of the number of objects in the sample for six different redshift un¬ 
certainties: Gz — 0.0001, 0.02, 0.04, 0.06, 0.08 and 0.1. The top 
row of plots show the fractional error for the correlation between 
environment and F-Low. The bottom row of plots show the frac¬ 
tional error for the correlation between environment and F-High. 
The left hand plots are for the iV* nearest neighbour method and 
the right hand plots are for the fixed aperture method. 

There are several points to note about the plots. Firstly the 
plots show that the fractional error decreases as the number of ob¬ 
jects in the sample increases. Secondly for a particular environment 
correlation (strong or weak) and and for a particular method, larger 
redshift uncertainties have larger fractional errors. Thirdly for a par¬ 
ticular correlation (strong or weak) the fractional error for the spec¬ 
troscopic redshift uncertainty as a function of the number of objects 
in the sample is approximately the same for the nearest neigh¬ 
bour and fixed aperture methods. Fourthly by comparing plots a) 
and b) for the weak correlation to plots c) and d) for the strong cor¬ 
relation it is clear that a larger spectroscopic trend correlation leads 
to a smaller fractional error for the same redshift uncertainty and 


this is true for both the A/^* nearest neighbour method and the fixed 
aperture method. 

We used the fractional error curves to obtain the number of 
spectroscopic objects required to achieve a particular fractional er¬ 
ror. In this way we obtained the ratio of the number of photometric 
objects to the number of spectroscopic objects - the sample size 
factor, i.e how many times larger the photometric sample must be 
relative to the spectroscopic sample to achieve the same fraction 
error. 

Figure [16] shows four plots of the sample size factor as a func¬ 
tion of the number of spectroscopic objects for five different red¬ 
shift uncertainties: Gz = 0.02, 0.04, 0.06, 0.08 and 0.1. The ar¬ 
rangement of the plots is the same as for Figure [15] The top row of 
plots are for the correlation between environment and F-Low. The 
bottom row of plots are for the correlation between environment 
and F-High. The left hand plots are for the A^* nearest neighbour 
method and the right hand plots are for the fixed aperture method. 
The sample size factor is fairly insensitive to the number of objects 
in the sample as the curves are nearly horizontal. The sample size 
factor increases with redshift uncertainty. 

At a redshift uncertainty of 0.1 approximately 6 — 16 times 
more photometric objects than spectroscopic objects are required 
to measure environment correlations with equivalent fractional er¬ 
rors for spectroscopic correlations of 0.3 —0.7 respectively. Galaxy 
properties that have a strong spectroscopic environment depen¬ 
dence are easier to detect in absolute terms with photometry than 
properties that correlate weakly with the spectroscopic environ¬ 
ment measurements. However we calculate the sample size factor 
required to obtain a photometric measurement of the correlation 
with the same fractional error as was obtained with spectroscopy. 
We find that galaxy properties that have a strong spectroscopic en¬ 
vironment dependence actually require a larger sample size factor 
to obtain the same fractional error than galaxy properties that de¬ 
pend weakly on spectroscopic environment. This is because spec¬ 
troscopic measurements of strong correlations have a small error 
whereas the uncertain redshifts still contribute a large error to pho¬ 
tometric measurements. Even more photometric objects are there¬ 
fore required to bring the photometric error down to a competitive 
level. 

We note that the sample size factors for the A^* nearest neigh¬ 
bour method and the fixed aperture method are very similar and 
so there is no preferred method. Excluding considerations of the 
dependence of galaxy properties on mass we therefore find that rel¬ 
atively small sample size factors are required to make equivalent 
detections of the correlations between galaxy properties and envi¬ 
ronment. 

5.3 Brief comparison with literature 

We complete this sectio n with a few paragra phs to compare this 
work with the literature. ICooner et al.l (l2005h presented the most 
extensive range of tests to date regarding measures of galaxy en¬ 
vironment and the impact of photometric redshifts. There are a 
numbe r of important differences between the work bv lCooner et al.l 
([2005h and this work. Using mock galaxy catalogues to mimic 
galaxy surveys Cooper et al. investigated a higher and larger red¬ 
shift range: 0.7 < z < 1.4 and considered a 1 squared degree field 
that contained ~22,000 galaxies. We have studied a low redshift 
range: 0.02 < z < 0.085 using data extracted from the SDSS DR7 
which co vers nearly 7000 squ ared degrees and includes ~ 150,000 
galaxies. ICooper et ahl (l2005h concluded that photometric redshift 
measurements with errors > 0.02 severely limit studies of the lo- 
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Figure 15. Shows the fractional error of the Spearman rank correlation coefficient between the environments and the faked data where the correlation of the 
faked data to the spectroscopic benchmark environments was constructed to be 0.3 (a and b) and 0.7 (c and d) for the nearest neighbour method (a and c) 
and the fixed aperture method (b and d) as a function of the number of objects in the sample for six redshift uncertainties: az = 0.0001, 0.02, 0.04, 0.06, 0.08 
and 0.1. 


cal galaxy environment. Nevertheless they found a Spearman rank 
correlation coefficient of ~0.3 between real space and photomet¬ 
ric measures_ofgatey environment for cr^ = 0.05 (see Table 3 
in ICooper et al.ll200^ . This is non-zero and so suggests that there 
actually is still a me asurable environmen t signal, albeit a weak 
one. We point out that ICooper et al.l (l2005h employed velocity cuts 
< 2000 km/s in their analysis and argued that deeper apertures 
w ould not measure the local environment. We build on the work 
of ICooper et al.l (l2005h by explicitly testing deeper apertures and 
use velocity cuts of 1000 — 20, 000 km/s. We show that the corre¬ 
lation between environments based on photometric redshifts and a 
set of spectroscopic benchmark environments depends on the aper¬ 
ture parameters. To obtain the optimal correlation we find that the 
aperture depth must increase with redshift uncertainty. We argue 
that using deeper apertures is a measurement integrated over the 
large scale structure (along the line of sight) however it is also a 
proxy for the local galaxy environment in photometric surveys. Re¬ 
ferring to Figure [T0|in this work we estimate the optimal Spearman 
rank correlation coefficie nt at cr;^ — 0.05 to b e -^0.5. This is larger 
than what was reported in ICooper et al.l (l2005h . However the differ¬ 


ence can be explained by referring to plots g) and h) in Figures 8 
or 9. Ado pting a non-optimal velocity cut (i.e. ~2000 km/s) as was 
chosen in ICooper et ahl (l2005h results in a degraded Spearman rank 
correlation coefficient. 

We also comp are our work with a more recent study by 
iFossati et^ (l2015h . This study is based on semi-analytic models 
of galaxy formation and inves tigates a redshift snapshot at z^l. 
Figure 4 in iFossati et al.l (l2015h shows the galaxy density measured 
within cylindrical apertures with radii of 0.75 Mpc for high-res 
measurements (i.e. simulating high quality spectroscopic redshift 
measurements) and for photometric redshift measurements (ob¬ 
tained by convolving the high-res measurements with a Gaussian 
distribution l eading to a redshift accuracy of ~0.015 at z~1.0). In 
the study bv iFossati et al.l (l2015h velocity cuts of 1500 km/s and 
7000 km/s were used for the high-res measurements and photo¬ 
metric redshift measurements respectively. Although there are dif¬ 
ferences in redshift and aperture size we compare their Figure 4 
with our Fi gure [5] In a qualita tive sense the plots are similar. In 
Figure 4 bv lFoss^i et al.l (l2015h the photometric measurements of 
the high densities recover the high-res density measurements as the 
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Figure 16. Shows the number of times larger than the spectroscopic sample the photometric sample must be to obtain equivalent fractional errors for the 
Spearman rank correlation coefficient between environment and faked data where the correlation of the faked data to the spectroscopic benchmark environments 
was constructed to be 0.3 (a and b) and 0.7 (c and d) for the nearest neighbour method (a and c) and the fixed aperture method (b and d) as a function of 
the number of spectroscopic objects for five redshift uncertainties: az = 0.02, 0.04, 0.06, 0.08 and 0.1. 


blue contours follow the dashed 1-to-l line. We also recover the 
high density environments better than the low density environments 
however our photometric environment measurements fall short of 
the 1-to-l line. This is because we c ompute a volume de nsity nor¬ 
malized using the mean density and iFossati etaP (12015h presents 
a surface density. Both studies report a reduced dynamic range for 
the photometric measurements compared to the spectroscopic mea¬ 
surements. 

iFossati etaP(l2015h also investigate trends with galaxy envi¬ 
ronment. They present a method to “clean” the population of low 
mass cen trals and clas s ify th e central and satellite galaxies. Fig¬ 
ure 10 in IFossati et al.l (l2015h shows the red fraction for the cen¬ 
tral galaxies as a function of the surface density and stellar mass 
for the high-res and photometric redshift measurements. The trend 
obtained with the high-res measurements is reproduced using the 
photometric measurements; alt hough this conclusion comes with 
caveats (^see IFossati et ahlEoisL for more details). In this paper we 
also report that we are able to reproduce trends, albeit some degra¬ 
dation, using photometric redshifts but instead of using simulations 
we use observed data from the SDSS. We present the red fraction as 


a function of mass and environment JPeng et al.ll201Qh using pho¬ 
tometric redshifts (see the second row of Figure [TJ]) 

Lastly we add to earlier works by attempting to give a “ball¬ 
park” estimate of how much larger than spectroscopic samples pho¬ 
tometric redshift samples must be to measure environment correla¬ 
tions. 


6 CONCLUSIONS 

The next generation of photometric surveys such as the DBS and 
Euclid will deliver enormous datasets. The cost of surveying enor¬ 
mous cosmological volumes is redshift precision. The aim of this 
paper is to examine the impact of uncertain redshifts on the mea¬ 
surement of galaxy environment. 

We use a low redshift (0.02 < z < 0.085) sample of galax¬ 
ies from the seventh data release of the SDSS and investigate two 
methods to measure galaxy environment: (i) the iV* nearest neigh¬ 
bour method (ii) and the fixed aperture method. We study a range 
of aperture parameters. 
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For the N'^ nearest neighbour we study = 1 —10. This cor- 
responds to median comoving projected scales of 0.3 — 3.4 Mpc. 
To ensure that the projected size of the apertures are well matched 
between the two methods we use fixed apertures with radii of 
0.1 — 3.0 Mpc centred on the galaxy targets. For both methods we 
employ a range of velocity cuts: 1000 — 20,000 km/s. A velocity 
cut of 20, 000 km/s extends through the entire redshift range. The 
possible comoving depths of the apertures are 27.9 — 1271.7 Mpc. 

We establish a set of spectroscopic benchmark environments 
for the iV* nearest neighbour method and the fixed aperture 
method. The aperture parameters for these benchmark measure¬ 
ments are A/ = 4 and dv = 1000 km/s for the iV* nearest neigh¬ 
bour method and r = 1.8 Mpc and dv = 1000 km/s for the fixed 
aperture method. 

We calculate environment measurements for a spectroscopic 
and photometric redshift sample and find the redshift uncertainty 
of the SDSS photometric redshift sample is 0.0185. We find that 
the A/* nearest neighbour method gives a smoother environment 
distribution compared to the fixed aperture method, especially for 
sparse environments. The A/* nearest neighbour method typically 
also has a larger dynamic range than the fixed aperture method. 

We compare the spectroscopic environment distributions to 
the photometric environment distributions and find for both meth¬ 
ods that the photometric range is smaller than the spectroscopic 
range for fixed aperture parameter values. We attribute the reduc¬ 
tion in dynamic range to the uncertain photometric redshifts. Un¬ 
certain redshift measurements tend to scatter galaxies away from 
dense regions and into less dense regions. 

We study the impact of the aperture parameter values by cal¬ 
culating the Spearman rank correlation coefficient between envi¬ 
ronment measurements and the spectroscopic benchmark environ¬ 
ments. 

We construct a set of simulated photometric redshift cata¬ 
logues by displacing the spectroscopic redshifts and correcting the 
r-band absolute magnitudes. The catalogues cover a range of red¬ 
shift uncertainties: az = 0.0025 — 0.1. We use the simulated 
photometric redshift catalogues to study the impact of uncertain 
redshifts on the Spearman rank correlation coefficient between the 
photometric environments and the spectroscopic benchmark envi¬ 
ronments. 

We find that the Spearman rank correlation coefficient be¬ 
tween the photometric environments and the spectroscopic bench¬ 
mark measurements generally decreases as the redshift uncertainty 
increases. The correlation depends also on the choice of aperture 
parameter values. For low redshift uncertainties {az < 0.02) there 
exists a relative ‘sweet spot’. For larger redshift uncertainties pro¬ 
viding the aperture parameter are large enough (e.g. A/ > 3, 
r > 1.2 Mpc, dv > 10,000 km/s and smaller than the maxi¬ 
mum values we study) the correlation is relatively insensitive to the 
aperture parameters. 

In addition we find the parameter values that result in the 
strongest correlation between the environments and the spectro¬ 
scopic benchmark measurements at each redshift uncertainty. The 
parameters that control the projected size of the apertures i.e. N 
for the A/* nearest neighbour method and r for the fixed aperture 
method are fairly insensitive to the redshift uncertainty. The op¬ 
timal values for these parameters are essentially the same as the 
benchmark values. The optimal value for the velocity cut however 
increases with redshift uncertainty. This is to ensure that galaxies 
that are increasingly scattered along the line of sight by the uncer¬ 
tain redshift measurements are capture within the apertures. 

Studying the Spearman rank correlation coefficient as a func¬ 


tion of the redshift uncertainty using the optimal parameter val¬ 
ues at each redshift uncertainty we find that the correlation decays 
rapidly at low uncertainties {az < 0.02) and then more gradually 
at larger redshift uncertainties {az > 0.02) 

Adopting the optimal parameter values at a redshift uncer¬ 
tainty of 0.1 we find that Spearman rank correlation coefficient 
between the photometric environments and the benchmark envi¬ 
ronments is ~0.4 for both the A/* nearest neighbour method and 
the fixed aperture method. 

As a case study we examine the bivariate dependence of the 
red fraction of galaxies on log(mass) and environment. Using the 
spectroscopic redshift measurements and the A/* nearest neighbour 
method {N = 4 and dv — 10 00 km/s) we reproduce the result 
presented by IPeng et aP (l2010h . In addition we present red frac¬ 
tion surfaces using the fixed aperture method (r = 1.8 Mpc and 
dv — 1000 km/s). We also study the red fraction bivariate de¬ 
pendence with the SDSS photometric redshifts and are able to ob¬ 
tain similar, albeit slightly deteriorated surfaces. The information 
to split the red fraction by mass and environment is still present in 
the SDSS photometric redshift measurements. Using the simulated 
photometric redshift catalogues and the optimal parameter values 
we find that the deterioration of the red fraction surfaces increases 
with redshift uncertainty. Nevertheless for the A/* nearest neigh¬ 
bour method even with large redshift uncertainties {az — 0.06) 
and including an additional scatter of 0.3 dex in the mass estimates 
the environmental dependence of the red fraction does not break 
down entirely. 

Using a Monte-Carlo type approach we calculate how much 
larger, than a spectroscopic sample, a photometric sample needs to 
be to obtain equivalent fractional errors for trends between galaxy 
properties and galaxy environment. To do this we use sets of faked 
data that correlate with the spectroscopic benchmark environments 
but not with mass. We define the fractional error as the ratio of the 
median Spearman rank correlation coefficient between the faked 
data and environment and the standard deviation of the correlation 
coefficient distribution from a set of 224 realizations. 

We find that photometric redshift samples with the expected 
DBS redshift uncertainty of 0.1 must be approximately 6 — 16 times 
larger than spectroscopic samples to measure environment correla¬ 
tions with equivalent fractional errors assuming that the correlation 
between the galaxy properties and the spectroscopic environment 
measurements was 0.3 — 0.7 respectively. The performance of the 

nearest neighbour and the fixed aperture method are similar 
when the optimal parameter values are adopted. 

This study was limited to a low redshift sample and the results 
may not hold at higher redshifts. Nevertheless the largest spectro¬ 
scopic pencil beam surveys, such as DEEP2 and COSMOS only 
cover a few squared degrees on the sky. DBS is expected to measure 
5000 square degrees and Euclid will be a full sky survey. Hence 
such sample size factors are easily achievable in the next generation 
of large scale photometric surveys potentially enabling competitive 
measurements of the correlations between galaxy properties and 
their environments. 
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